Signatures of orbital loop currents in the spatially resolved local density of states 
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Polarized neutron scattering measurements have suggested that intra-unit cell antiferromagnetism 
may be associated with the pseudogap phase. Assuming that loop current order is responsible for the 
observed magnetism, we calculate some signatures of such circulating currents in the local density 
of states around a single non-magnetic impurity in a coexistence phase with superconductivity. We 
find a distinct C4 symmetry breaking near the disorder which is also detectable in the resulting 
quasi-particle interference patterns. 
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I. INTRODUCTION 

The understanding of the pseudogap phase in cuprate 
materials constitutes an outstanding challenge in con- 
densed matter physics. At present, it remains unsolved 
whether the phase is caused by singlet formation natu- 
rally exhibited in strong-coupling models, or rather by 
an elusive spontaneously symmetry-broken phase. 

In a recent series of spin-polarized neutron scat- 
tering experiments, it has been observed that intra- 
unit cell antiferromagnetic order sets in at T*, the 
pseudogap temperature in the underdoped regime of 
YBa2Cu306+x and HgBa2Cu04+5.^^ Recently, a short- 
range magnetic intra-unit cell signal was also identified 
in Lai.9i5Sro.o85Cu04 though at temperatures signifi- 
cantly below T*.^ These results suggest that the pseu- 
dogap phase may in fact be associated with a true phase 
transition with an associated spontaneously broken sym- 
metry. Varma has predicted that the pseudogap phase 
is caused by equilibrium loop currents which break the 
time-reversal symmetry but preserve the translational 
invariance of the Cu02 lattice. ^'^ A recent theoretical 
study^ found that such loop currents can be stabilized 
by using the apical oxygen ions leading to ordered mo- 
ments rotated away from the axis perpendicular to the 
Cu02 plane in qualitative agreement with the neutron 
measurements. 

The possible existence of loop order remains highly 
controversial, however, since local probes fail to detect 
the weak magnetic field associated with such orbital cur- 
rents. For example, nuclear magnetic resonance (NMR), 
nuclear quadropole resonance (NQR), and muon spin re- 
laxation (/iSR) have failed to detect any onset of mag- 
netic order setting in at T*.^^^ Most recently, a high- 
precision zero-field /iSR study on La2-xSra^Cu04 found 
no evidence of static magnetic order in the pseudogap 
phase at doping levels above x = 0.13, which is outside 
the spin-glass phase in LSCO.^^ Various proposals, in- 
cluding minority phases, loop-order fluctuations, and lo- 
cal muon destruction of loop order have surfaced to rec- 
oncile the apparent contradiction between local probes 
and neutron experiments. ^^'^^ At present, the possible 



existence of loop order remains unsettled and therefore 
new proposals for its detection seem desirable. 

Scanning tunneling spectroscopy (STS) is a natural ap- 
proach to take. STS experiments on Bi2Sr2CaCu208, 
for example, have found evidence for intra-unit cell "ne- 
matic" symmetry breaking, which shows up as an in- 
equivalency of oxygen p^^ and p^ orbitals.^^ Unfortu- 
nately, loop order does not have such a simple experi- 
mental signature, since the local density of states (LDOS) 
has the same C4 symmetry as the lattice. The loop cur- 
rents are expected to generate a nonzero density of states 
at the Fermi energy, ^^ but are not unique in this; other 
mechanisms, such as disorder, also generate low energy 
excitations. 

In this paper, we show that there is a clear signature 
of loop currents if one considers the LDOS near an iso- 
lated impurity. In the literature, disorder effects have 
been previously suggested as a probe of the so-called 
d-density wave candidate for the pseudogap phase. ^^^^^ 
Here, our aim is not to answer whether loop currents ex- 
ist, but rather to suggest new "smoking gun" experiments 
based on their assumed existence. We focus on the low- 
temperature coexistence phase of loop order and d-wave 
superconductivity since this is where STS measurements 
are most readily performed. We use a mean-field three- 
band model to study the effects of loop order on the 
LDOS near non-magnetic impurities. Our main result 
[Fig. 1] is that there is a distinct C4 symmetry breaking 
in the vicinity of the impurity, and an associated splitting 
of the so-called "octet" vectors in the Fourier transformed 
scanning tunneling spectroscopy (FT-STS) patterns aris- 
ing from quasi-particle interference (QPI). These effects 
can be explained solely from symmetry arguments, and 
result from momentum-selective shifts of the nodal Dirac 
points in the d-wave superconductor in the presence of 
loop order. 



II. MODEL 

We employ the three-band Hubbard Hamiltonian with 
mean- field decoupled interactions.^ As shown recently. 



loop currents are stable in this model when the nearest- 
neighbor repulsion Vpd is sufficiently large compared 
to the onsite repulsions on the copper Ud and oxy- 
gen Up sites. ^^ In this paper, we assume the exis- 
tence of loop currents and study their detectable con- 
sequences. The three-band Hamiltonian is a minimal 
model containing the essential symmetries of the prob- 
lem. Experiments^^ suggest that the current loops in- 
volve out-of-plane orbitals,^ which are not part of the 
model. These additional orbitals are neglected for sim- 
plicity, as they do not change the rotational symmetries 
of the Hamiltonian, and it is these symmetries that de- 
termine the symmetry of the LDOS pattern near an im- 
purity. 

In the three-orbital basis '^^(k) 



(4'i^l,k'4,k) 



the normal-state Hamiltonian including orbital currents 
takes the form 



n = ^^H^){ H'{k)^^H^{k)) ^{k), (1) 

where H^{k) and H^{k) are 3 x 3-matrices. H^{k) de- 
scribes the hopping part of the Hamiltonian as well as 
the decoupled on-site interactions, whereas H^(k) car- 
ries the decoupled oxygen-copper interactions and is thus 
responsible for the loop currents. 

Writing out equation (1) yields the following expres- 
sion for the Hamiltonian: 
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where q = cos{ki/2) and Si = sm{ki/2). Following e.g. 
Hybertsen et al.^^^ we use the following set of parameter 
values; tpd = -1.3, V = OMpd, e^ = -1.5, Cp = -5.0 
(all in eV). In Eq.(2), R denotes the mean- field order pa- 
rameter for the orbital currents, and is set to either zero 
or, in the current-carrying case, to 0.1 i This phase choice 
of R leads to a current-carrying state in accordance with 
the original mean- field formulation by Varma.^ The main 
results of this paper are not sensitive to the amplitude of 

To include d-wave superconductivity, we change 
to a six-operator Nambu basis given by '^(k) = 

(^kt'^l,kt'^y,kt'^-H':Pa^ -ki^P?/ -ki). The Hamiltonian 
is then given by 



n = J2^H^) 



F(k) A(k) 
At(k) -H*{-k) 



^(k), (3) 



where A(k) is the 3x3 matrix describing the Cooper- 
pairing in our model. For simplicity, we take all Cooper 
pairing to be on the copper orbital which in the present 
case is a good approximation because the band crossing 
the Fermi level is mainly of Cu orbital character when 
Cd = —1.5, Cp = —5.0. We assume the gap function to 
have (ia,2_^2 -symmetry with 

A(k) = Ao(cos(A:^) - cos(A:^)) • diag(l, 0, 0). 



(4) 



Other choices for A(k), involving pairing in the oxygen 
orbitals, are possible; however, these will not change the 
symmetry of the LDOS pattern. To allow us to focus on 
spectral features below the gap, we take an artificially 
large Aq = 0.25 eV. 



The free Green's function for this system is now, for 
each k, a 6 X 6 matrix in orbital and spin indices given 

by 



GO(k, w) = (l(w + irj) - H'{k))-\ 



(5) 



where i^'(k) is the 6x6 matrix of Eq.(3), and rj is an 
infinitesimal regulator. 

Finally, we introduce a pointlike non-magnetic impu- 
rity at the copper site belonging to the unit cell at tq. 
The impurity Hamiltonian is 

Himp = l^mp-(diag(l,0,0)©diag(-l,0,0))(5(r-ro). (6) 

Using a standard T-matrix formalism, the full Green's 
function is then given by 

G{t,t\uj) = G^(r-r^cJ)+G^(r-ro,cJ)THG°(ro-r^cJ), 

(7) 
where 



T{lo) = (1 - Hi^pG°{0,u)r^Hi^p. 



(8) 



From the full Green's function, the local density of states 
p^(r,cj) on the orbital i in the unit cell r is readily ob- 
tained from the formula (for the superconducting state) 



piir.uj) = -Im Y^ G^^,^^(r,r,acj)/7r, 



(9) 



a=± 



with a labeling spin. The QPI spectrum is then obtained 
from the Fourier transformed density of states: 

(10) 



III. RESULTS 

The LDOS patterns in the vicinity of the impurity are 
shown in Fig. 1, along with the corresponding current 
patterns. The current pattern denoted by ^ = 1 is gen- 
erated by the Hamiltonian (2). The other two current 
patterns are obtained by changing the off-diagonal en- 
tries of i^(k): 

A = 2: Hi20^)=2itpdS:,-Rc^, (11) 

His{k) = -2itpdSy^Rcy. (12) 

^ = 3 : Hi2{k) = 2itpdSx + Rsx. (13) 

Hi^{\^) = -2itpdSy^Rsy. (14) 

Schematically, the unit cell in real space can be repre- 
sented as 



Py 
d Px 



which also provides a legend for reading Fig. 1. As can be 
seen from Fig. 1, the symmetries of the LDOS patterns 
are inherited from the spatial symmetries of the loop or- 
der, with the obvious additional rule that the directions 
of the current arrows play no role in the LDOS. These 
symmetries are tabulated in Table I, and are easily ob- 
tained from the symmetries of Eq. (7). We stress that 
in the clean phase, the LDOS has the full symmetry of 
the lattice, and that it is only near an impurity that the 
distinctive signatures of the loop currents are revealed. 

We can understand the broken spatial symmetries pre- 
sented in Fig. 1 by consulting the contours of constant 
energy (CCE) of the bare band structure. As in the case 
of a pure d-wave superconductor, the band structure in- 
side the gap yields closed banana-shaped CCE which in- 
crease in size with the energy u as shown in Fig. 2(a). 
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FIG. 1: C4 symmetry breaking near an impurity positioned 
at the origin (0,0) by loop currents, (a)-(c) Possible loop 
current patterns involving a central Cu and four neighboring 
O orbitals. (d)-(f) Corresponding LDOS at cj = 0.2 eV on 
the Cu and O orbitals near a single nonmagnetic impurity at 
the central Cu site marked by a white cross. Black sites are 
vacant, numbers label Cu atoms. Results are for (a),(d) 
^=1, (b),(e)^ = 2, (c),(f)^ = 3. 
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No 
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Yes 


Yes 


Yes 
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Yes 


Yes 


Yes 


C4 


No 


No 


Yes 



TABLE L Symmetries of the LDOS near an impurity. 



In Fig. 2(b) we see the effect of orbital currents with 
^ = 1 on the CCE; while two of the Dirac cones are left 
invariant, the other two are shifted upwards/downwards 
in energy. ^^ These distorted bananas respect the same 
symmetries as the current patterns. Changing the cur- 
rent direction corresponds to interchanging the large 
and small bananas whereas changing the current pattern 
A = 1 ^ A = 2 causes a C4 rotation of the LDOS as 
expected. 

Experimentally, the CCE may be inferred from QPI 
patterns measured by tunneling experiments. For pure 
d-wave superconductors, these QPI patterns have been 
explained using a simple "octet" model, ^^^^ in which 
peaks in the QPI pattern are attributed to scattering be- 
tween portions of the CCE for which the joint density 
of states is high. Figure 2(a) shows the octet vector q^, 
which labels the most prominent "intra-banana" quasi- 
particle scattering process. 

The effect of ^ = 1 loop currents on q^ is illustrated 
in Fig. 2(b): the q^ peaks along the left diagonal are 
split into q^ and q^, while those along the right diagonal 
remain unchanged. This uniaxial distortion is certainly 
seen in the calculated QPI spectrum shown in Fig. 3(a- 
d) which display |p(q, cj)| at cj = 0.09 eV as a function 
of q. In addition. Fig. 3(e) shows how the identification 
of the q^ peaks in the QPI spectrum is in very good 
agreement with the octet model prediction obtained di- 
rectly from the CCEs. As is evident from e.g. Fig. 2(b), 
the other (than q^) characteristic q-vectors of the octet 
model connecting the banana tips will be slightly rotated 
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FIG. 2: Contours of constant energy for energies near 
cj = 0.2 eV (a) without and (b) with A—1 orbital loop 
order. 
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FIG. 3: Plot of |p(q,cj)| at cj = 0.09 eV (a) without and (b) 
with A—\ loop currents. Arrows indicate (a) q^ and (b) q^ 
and q^ peaks. Cuts of |p(q, cj)! are shown along qy — —Qx 
(solid) and qy = qx (dashed) (c) without and (d) with loop 
currents. Arrows indicate the same peaks as in (a) and (b). 
(e) Dispersion of the q^ peaks as a function of uj. There is 
good agreement between the peak positions in p(q, cj) 
(points) and the octet model (solid lines). 



off the center axis in the presence of circulating current 
loop order. 



IV. DISCUSSION AND CONCLUSION 



There are a number of possible caveats to the present 
proposal that should be discussed. First of all, the pres- 
ence of dynamical loop order fluctuations will clearly be 
detrimental to the observation of C4 spatial symmetry 
breaking. Second, the existence of various domains with 
respectively A = 1 and A = 2 loop order will "sym- 
metrize" the QPI images and roughly result in a smeared 
version of the QPI from a pure d-wave superconductor 
without current loops. The local C4 symmetry breaking 
near a single impurity should still be detectable in this 
case. For BSCCO, however, it is expected that the LDOS 
near single-site disorder is further complicated by intrin- 
sic disorder and the structural supermodulation.^^'^^ Fi- 
nally, we should also mention that the current calculation 
is not self-consistent and the feedback effect of the impu- 
rity on the local current loops is not included. The effect 
of a charged impurity on the loops has been addressed 
previously in the discussion of the potential disturbance 
of the /i+ on the current order. ^^'^^ In the present case 
it is unlikely that the local suppression of loop currents 
would restore the symmetry. Thus, we expect the overall 
conclusion of C4 symmetry breaking near disorder sites 
to remain valid. 

In summary, we have shown how spatial C4 symmetry 
may become broken in the LDOS around nonmagnetic 
impurities in the coexistence phase of loop order and d- 
wave superconductivity. In addition there exists an as- 
sociated splitting and/or rotation of the so-called octet 
peaks evident in the quasiparticle interference patterns. 
These predictions should be testable by future tunneling 
measurements. 
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